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ABSTRACT 

Using CHANDRA, we investigate the spatial temperature distribution of the intracluster medium (ICM) 
within 700 kpc of the center of the massive merging cluster MACSJ0717. 5+3745 at z = 0.55. Combining 
the X-ray evidence with information about the distribution and velocities of the cluster galaxies near the core 
provides us with a snapshot of the three-dimensional geometry and dynamics of one of the most complex 
cluster studied to date. We find MACSJ0717. 5+3745 to be an active triple merger with ICM temperatures ex- 
ceeding 20 keV. Although radial velocity information and X-ray/optical offsets indicate that all three mergers 
proceed along distinctly different directions, the partial alignment of the merger axes points to a common origin 
in the large-scale filament south-east of the cluster core. Clear decrements in the ICM temperature observed 
near two of these subclusters identify the respective X-ray surface brightness peaks as remnants of cool cores; 
the compactness and low temperature of 5.7 keV of one of these features suggest that the respective merger, a 
high-velocity collision at 3,000 km s , is still in its very early stages. Looking beyond the triple merger, we 
find the large-scale filament to not only provide a spatial as well as temporal arrow for the interpretation of the 
dynamics of the merger events near the cluster core; we also find tantalizing, if circumstantial, evidence for 
direct, large-scale heating of the ICM by contiguous infall of low-density gas from the filament. 
Subject headings: galaxies: clusters: individual (MACSJ07 17.5+3745); X-rays: galaxies: clusters 



1. INTRODUCTION 

Major cluster mergers are the most energetic events in the 
universe, and one of the most important phenomena for the 
formation of large-scale structure in the hierarchical scenario. 
The enormous energy released in cluster mergers accelerates 
particles to relativistic energies and can, in the most massive 
systems, heat the intracluster medi um (ICM) to temperatures 
exceeding 20 keV (see review bv lSarazinj|2002l) . In an astro- 
physical context, dynamically simple two-component merg- 
ers like 1ES0657-56 have been instrumental in obtaining the 
currently best c onstraint on the self-in t eraction cross sectio n 
of dark matter dMarkevitch et al.ll2004t iRandall et al] l2008)'. 
Few mergers, however, are as straightforward to interpret as 
1ES0657-56 which consists of only two well separated com- 
ponents, and whose merger axis is almost perpendicular to our 
line of sight. While the majority of massive mergers are thus 
among the most challenging systems for studies of the dynam- 
ics and interactions of gas, galaxies, and dark matter, they also 
offer the greatest opportunity to study with good statistics the 
full complexity of cluster mergers. 

A prime example of a complex, major cluster merger is 
MACSJ0717.5+3745 (z = 0.54 5). Discovered during th e 
Massive Cluster Survey (MACS lEbeling et ail 1200 lL 120071) . 
the system is one of the most X-ray luminous and best- 
studied massive clusters at intermediate redshifts. The wealth 
of data accumulated for it from X-ray to mm wavelengths 
make MACSJ0717. 5+3745 a promising target for a com- 
prehensive study of the full range of physical mechanisms 
at work during the assembly of massive clusters. At radio 
frequencies, synchroton emission from relativistic electrons 
in MACSJ07 17.5+3745 was detected in the form of both a 
radio halo and a very extended, steep-spectrum radio relic 
(lEdgeetal.1 120031 The same relativistic electrons are ex- 



1 see also Mahdavi et al. (2007); Bradac et al. 12008) for results obtained 
for A 520 and MACSJ0025.4-1222. 



pected to also produce hard X-ray photons via inverse Comp- 
ton scatterin g of Cosmic Microwav e Background photons, as 
reported by Petrosia n et all (l2006h . The spatial distribution 
of the galaxy po pulation in and aroun d MACSJ07 17.5+3745 
was studied by Ebeling et al] d2004h who found the clus- 
ter to be connected to an extended linear filamentary struc- 
ture, consistent with the picture seen in numerical simu- 
lations (and observa tionally supported by the findings of 
iKartaltepe et alj|2008l) that the most massive halos are more 
likel y to be surrounded by pro nounced large-scale structure 
(e.g. lColberg et alJ20 00. 1999). Finally, a spectroscopic anal- 
ysis by Ma et al. (2008) found strong evidence for star forma- 
tion in many galaxies in and around this cluster, likely caused 
by the ongoing growth of this complex system. In this let- 
ter, we study variations in the temperature of the ICM of 
MACSJ07 17.5+3745 and their correlation with features in the 
spatial and radial-velocity distribution of the cluster galaxy 
distribution, in order to assess the dynamical state and history 
of the cluster core. 

2. X-RAY DATA 

MACSJ07 17.5+3745 was observed with the Chandra Ad- 
vanced CCD Imagi ng Spectrometer Imaging array (ACIS-I; 
Garmire et al. 2003) in January 2003 for a total exposure time 
of 60 ks in VFaint mode (ObsID 4200). We applied standard 
Chandra data reduction procedures 2 to reprocess the event file 
using CIAO 3.3.0.1 and CALDB 3.2.1. Our spectroscopic 
analysis made use of a "blank-sky" background file 3 to re- 
move low-level soft X-ray emission covering a majority of 
the ACIS-I array. 

An adaptively smoothed representation of the data within 
the full field of view of ACIS-I is shown in Fig. Q] The 
disturbed X-ray morphology around the core of the cluster, 
which contains at least four X-ray surface-brightness peaks, 

2 |http : //cxc . harvard, edu/ciao /guides/ acis-data . htmlj 
3 see http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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FIG. 1 . — X-ray emission in the field of MACSJ07 17.5+3745 as observed 
with Chandra's ACIS-I detector in the 0.5-7 keV band. The image has been 
weighted by the exposure map and adaptively smoothed using the asmooth 
algorithm (Ebeling et al. 2006) requirin g a minimal significance of 99% wit h 
respect to the local background (see also Ma et al. 2008; Ebeling et al. 2007). 
The shaded region marks the chip gaps. Contours are logarithmically spaced 
by factors of two, starting at six times the value of the global X-ray back- 
ground. 

is immediately obvious in the upper right corner. In addi- 
tion, weak but detectable X-ray emission extends to the south- 
eastern corner of Fig.Q] matching th e large-scale filament dis- 
covered in the galaxy distribution by Ebelin get al.l d2004l) . 4 

3. OPTICAL DATA 

MACS J07 17.5+3745 was observed with the Advanced 
Camera for Surveys (ACS) aboard the Hubble Space Tele- 
scope (HST) in the F555W and F814W passbands on April 
2, 2004 for an effective exposure time of 4470 and 4560s, re- 
spectively, as part of a comprehensive HST survey of the most 
distant MACS clusters (GO-09722). Groundbased observa- 
tions of the galaxy population of the cluster were conducted 
with Keck-II/DEIMOS from 2004 to 2008 with the aim of 
probing large-scale cluster dynamics and examining the effect 
of environment on galaxy evolution (Ma et al. 2008). 

4. X-RAY TEMPERATURE MAP 

As is apparent from its complex X-ray morphology, 
MACSJ0717. 5+3745 is far from virialized, and the gas tem- 
perature can be expected to vary significantly in its core 5 . In 
order to derive a spatial temperature map, we need to de- 
fine suitable regions for X-ray spectroscopy 6 . Ideally, these 
regions should follow the temperature structure. The lat- 
ter, however, is a priori unknown, and the limited photon 
statistics of our data do not allow us to extract spectra in ar- 
bitrarily small bins. We address this problem by adopting 
a two-level binning method. As a first step, we subdivide 
the image into regions whose shape is guided by the X-ray 
surface-brightness distribution, and each of which contains 

4 This includes the satellite cluster 2 Mpc south-east of the main system; 
we measure its ICM temperature to be 2.7 ± 0.3 keV. 

5 A summary of the global properti es of the hot gas in 
MACSJ07 17.5+3745 is given bv lMa et alJ<2008T) . 

6 We perform all s pectral fitting with Sherpa using a MEKAL plasma 
model I Mewe et al. 1985), with the absorption term fixed at the Galactic value 
of nH = 7.11 X 10 20 ~cm~ 2 , the redshift of the cluster set to z ci = 0.5446 
)Ma et aU2008T) . and abundance Z = O.3Z . 



TABLE 1 
Results from spectral fitting 



Region 


kT* 


x 2 


Binning Table 




(keV) 






B 1 


24.3 ± 5.3 


0.78 


Al, A5, A7 


B2 


16.1 ±4.0 


0.57 


A2, A3, A4, A8 


B3 


14.8 ± 2.9 


0.56 


A6, A9, A15 


B4 


12.9 ±4.1 


0.69 


A10,A11, A13 


B5 


10.9 ± 1.7 


0.71 


A20, A27, A28 


B6 


10.2 ± 2.5 


0.73 


A12,A17, A19 


B7 


10.2 ±2.4 


0.71 


A14, A18, A26 


B8 


8.7 ± 1.2 


0.84 


A16, A23, A32 


B9 


8.4 ±3.6 


0.70 


A22 


B10 


8.2 ± 1.5 


0.69 


A21, A24, A25 


Bll 


7.3 ±0.7 


0.73 


A29, A30, A31, A34 


B12 


5.7 ±0.8 


0.93 


A33 


The quote( 


errors denote la uncertainty. 



enough photons to allow a crude temperature measurement. 
We then merge adjacent regions with similar temperature in 
order to obtain a clearer view of temperature variations on 
larger scales. This procedure ensures that the final regions 
represent the temperature structure fairly, and also that the re- 
sulting measurements are statistically robust. 

4.1. Low-fidelity temperature map 

The initial regions for this process are determined b y run- 
ning the "accumulative binning" algorithm contbin (Sanders 
2006) on the exposure-weighted image 7 . We chose contbin 
because the algorithm's choice of bins is physically moti- 
vated, in as much as contbin attempts to follow the X-ray 
surface brightness, which is related to the ICM temperature 
and density. If the initial regions are sufficiently small, their 
choice will, by the end of the merging process, not affect the 
final temperature map. 

The results of this first-level binning step are shown in the 
top two panels of Fig. [2] As expected, the relatively low SNR 
chosen leads to large errors in the measured temperatures and 
to low x 2 values of ~ 0.6. Although the best-fit tempera- 
tures for some regions, mostly the hottest ones, are poorly 
constrained, we already see clear variations and trends of the 
ICM temperature within our study region. 

4.2. Re-binning 

As a second step, we examined the temperature map from 
Stage 1 and combined adjacent regions with similar temper- 
ature in order to obtain a clearer view of temperature varia- 
tions on larger scales, the goal being to reach at least 1200 net 
counts per bin by combining three (occasionally four) of the 
initial regions. The results of this rebinning process are shown 
in the bottom two panels of Fig. |2j the individual regions and 
their temperatures are listed in Table Q] 

Note that two of the initial regions were excluded from 
the rebinning and left unchanged, namely the ones coincid- 
ing with the main peaks in the X-ray surface brightness dis- 
tribution, A22 (B9) and A33 (B12). Both of these exhibit 
markedly cooler temperatures than their surroundings, mak- 
ing them likely candidates for cool cores of clusters still in 
the process of merging with the main system. This interpreta- 
tion is particularly plausible for the region around the global 

7 In contbin the image was smoothed to a signal-to-nois e ratio (SNR) of 
10, and we set the "constraint parameter" C i Sanders 2006) to 1.5 to prevent 
the bins from becoming too elongated. The SNR required of the resulting re- 
gions was 20, which corresponds to 400 net counts in regions with negligible 
background. 
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FIG. 2. — Left panels: X-ray temperature maps from our two-step binning process. The temperature is coded as indicated by the color bar. Right panels: 
the best-fit temperature and reduced x 2 value for each bin. Top: results for the initial regions generated by contbin. Bottom: results for regions created in the 
subsequent rebinning process. The second temperature shown for region Bl (open circle) is the result obtained when the embedded cool region B9 is included. 



peak in the X-ray emission from MACSJ0717. 5+3745, A33 
(B12), which is found to contain the coolest gas anywhere 
within our study region (kT = 5.7 ± 0.8 keV). By compar- 
ison, the temperature measured for the second putative cool 
core, region A22 (B9), is much higher and less well con- 
strained (kT = 8.4 ± 3.6 keV), but still in stark contrast to 
the much higher temperatures above 15 keV measured in its 
immediate vicinity (regions Al, A5, and A7). With all of 
these regions being viewed in projection, the true temperature 
difference is likely to be yet higher. Better constraints from 
a two-phase plasma model, especially relevant for region B9, 
would require significantly better photon statistics. 

Although the selection of regions to merge is to some de- 
gree subjective, we convinced ourselves that our final results 
and conclusions are robust. For instance, our perhaps most de- 
batable choice, the exclusion of the aforementioned cool core 
A22 (B9) from the merging of the surrounding regions Al, 
A5, and A7 (Bl), predictably increases the temperature mea- 
sured for Bl (kT = 24.3 ± 5.3 keV); including it, however, 
still leads to a statistically consistent result of kT = 21.0±5.4 
keV(Fig.|3. 

5. CLUSTER GALAXY DISTRIBUTION AND DYNAMICS 



Our tentative interpretation of the physical origins of the 
complex temperature and gas density distribution presented 
in the previous section can be tested and extended by study- 
ing the distribution and dynamics of the cluster galaxies in 
the same region. Shown in Fig. [3] is an overlay of the X-ray 
contours from Fig. |2]onto the HST image of the cluster core. 
We identify four major concentrations in the integrated opti- 
cal light, which is dominated by the large cluster ellipticals. 
Average redshifts and velocity dispersions for galaxies within 
the four circles marked in Fig.[3]are listed in Table[2] Unlike 
their counterparts in regions A, C, and D, the galaxies associ- 
ated with feature B are found to move through the cluster core 
at a very high relative radial velocity of over 3,000 km s , 
close to the maximal velocity expected from infall from infin- 
ity. No significant differences are found between the velocity 
dispersions measured for the four regions, although the results 
suggest that the structure in region C is the most massive one. 

6. INTERPRETATION 

A comparison of the gas and galaxy distribution shown in 
Fig.[3]as well as of the radial velocity information summarized 
in Table |2] yields important insights into the dynamics of the 
multiple components of the cluster core. Specifically, we note 
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FIG. 3. — X-ray contours as in Fig.[2]overlaid on the F555W/F814W optical 
image of the cluster core obtained with HST/ACS. Yellow circles centered on 
the four main peaks in the cluster light distribution mark the regions within 
which average radial velocities and velocity dispersions were computed for 
the enclosed galaxies (see Table|2j- 



TABLE 2 

Average radial velocities and velocity 
dispersions for galaxies in the regions 

MARKED IN FlG.[3] 



Label 


Num of redshifts 


AV* 
(km/sec) 


&' 

(km/sec) 


A 
B 
C 
D 


10 
9 
10 
7 


r, 7s +295 
z '°-339 

OQQO+252 

OiOo 242 

oni+843 

OJ± -8nn 


1761+ 234 



. _ : """-snr 

Calculated using the bootstrap estimator in ] 



iBeers et al] 



that, of the four galaxy concentrations discussed in Section|5J 
only A and C are reasonably well aligned with peaks in the X- 
ray surface brightness. The good X-ray/optical alignment of 
component C, combined with a systemic velocity and velocity 
dispersion consistent with those of MACS J07 17. 5+3745 as a 
whole, makes this subsystem a likely candidate for the actual, 
if disturbed, core of the main cluster. 

The dominant galaxies in regions B and D, on the other 
hand, are offset by about 25 arcsec (~160 kpc at the clus- 
ter redshift) from the respective X-ray peak - in both cases 
in a direction that points toward the filament to the SE 
of the main cluster. Offsets of this kind between gas 
and galaxies (and dark matter) have been observed be- 
fore in merging clusters, examples be ing the Bullet Clus- 
ter 1E S0657-56 dMarkevitch et al.l2004l) andMACSJ0025.4- 
1222 (IBradac et al. 1 120081) . and are readily explained by the 
difference in the self-interaction cross sections of the various 
cluster constituents. The X-ray temperature maps of Fig. |2] 



show the gas near regions B and D being cool compared to 
their surroundings (B9 and B12); we can thus safely identify 
them as the cool cores of clusters merging with the main body 
of MACSJ0717.5+3745 and falling behind as their respective 
galaxy populations advance. While the merger axis of system 
D is consistent with lying in the plane of the sky (and aligned 
with the large-scale filament prominently visible on the left 
of Fig. [T]), the combination of the observed X-ray /optical off- 
set with a very high radial velocity indicates that system B is 
moving through the cluster core along an axis that is much 
more inclined toward our line of sight. We note also that the 
high-velocity collision of component B most likely happened 
only recently, given that its cool core (region B12 in Fig. O 
appears to be still largely intact. 

A similar argument can be made to explain the observa- 
tional evidence in region A. The X-ray/optical offset suggests 
infall from the NW, i.e. in the direction of the filament. The 
trailing peak in the gas density is highly elongated though, and 
no cool gas is detected in the area, all of which suggests that 
we are witnessing back-infall of a cluster that originated from 
the filament and has already passed through the main cluster 
once. 

While the combined X-ray/optical evidence allows a 
straightforward and self-consistent interpretation of the dy- 
namics and three-dimensional geometry of the three seperate 
ongoing merger events described above, the gas distribution 
and temperature structure near the cluster-filament interface 
- the south-eastern corner of Figs. [2] and [3]- defies simple 
explanation. Although the extremely high temperatures mea- 
sured in regions Bl and B2 could to some extent be caused 
by shock heating from the ongoing mergers of components B 
and D, this explanation appears implausible for the entirety of 
the hot interface region extending from B3 to B4. The sheer 
size of this feature (about 1 Mpc), its alignment perpendicu- 
lar to the cluster-filament axis, and the simultaneous increase 
in gas temperature and X-ray surface brightness all point to 
a direct large-scale interaction between the filament and the 
cluster core. We therefore speculate that the prominent jump 
in gas density and temperature along the cluster-filament in- 
terface is the result of contiguous accretion of gas along the 
large-scale filament feeding MACSJ0717.5+3745. A deeper 
Chandra observation would allow this hypothesis to be tested 
and may yield the first credible characterization of the gas 
content of large-scale filaments. 
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